PROBING THE FLYBY ANOMALY WITH THE GALILEO CONSTELLATION 
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ABSTRACT 



In the last few years, the so-called flyby anomaly has been 
widely discussed, but remains still an illusive topic. This 
is due to the harsh conditions experienced during an Earth 
flyby as well as due to the limited data available. In this 
work, we assess the possibility of confirming and charac- 
terizing this anomaly by resorting to the scientific capa- 
bilities of the future Galileo constellation. 
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1. INTRODUCTION - THE FLYBY ANOMALY 



During the past couple of decades, a few deep-space 
probes that used an Earth flyby have apparently displayed 
an unexpected velocity change after their gravitational as- 
sist. This has become known as the flyby anomaly. 

This variation in the velocity was detected in the residu- 
als of the analysis performed to the Doppler and ranging 
data — which showed the impossibility of fitting the tra- 
jectory with a single hyperbolic arc, but allowed for a 
separate fit of the inward and outward paths, if an unex- 
pected velocity shift is introduced. It is highly localized 
at the perigee, where tracking through the Deep Space 
Network (DSN) is not available (with an approximate 4 h 
gap). The spatial resolution of the available reconstruc- 
tions, resulting form the 10 s interval tracking, does not 
allow for an accurate characterization of the effect, so that 
no corresponding acceleration profile exists; thus, only 
the variation of the probes' velocity (vis-a-vis kinetic en- 
ergy) is known. 

This flyby anomaly has so far been observed in Galileo, 



NEAR, Rosetta, and Cassini missions |[T]. A summary of 
the Earth flybys observed since the 1990's is shown in Ta- 
ble [1] A detailed discussion of the two Galileo (1990 and 
1992) and the NEAR (1998) gravity assists is presented in 
Ref. [2 |, including an analysis of the three earliest flybys 
where the anomaly was observed; it included an estimate 
of the accelerations generated by different known effects, 
in an attempt to single out possible error sources. 

It is estimated that the average acceleration associated 
to the flyby anomaly is of the order of 10^"* m/s^ lU. 
This is measured against the Earth oblateness, other solar 
system bodies, relativistic corrections, atmospheric drag. 
Earth albedo and infrared emissions, ocean tides, solar 
pressure, etc. . 

Subsequently, Ref. (J] extended this discussion to other 
possible error sources, comparing this 10^"* m/s^ fig- 
ure with several possible origins for additional unac- 
counted for accelerations: these include the atmosphere, 
ocean tides, solid tides, spacecraft charging, magnetic 
moments. Earth albedo, solar wind and spin-rotation cou- 
pling. The authors conclude that all of the considered 
effects are several orders of magnitude below the flyby 
anomaly. 

From these references, one can compile Table |2] for a 
quick overview of all effects. Taking into account a value 
of 10~^ m/s^ for the hypothetical flyby anomaly, all 
listed effects, except the Earth oblateness have lower or- 
ders of magnitude. This raises the issue of possible errors 
in the gravitational model of Earth: however, attempts to 
solve the flyby problem by changing the related second 
dynamic form factor J2 have yielded unreasonable solu- 
tions, and are unable to account for all flybys [[ZJ. 

An empirical formula to fit the flyby parameters has been 
proposed as a function of the declinations of the incom- 
ing and outgoing asymptotic velocity vectors, 6i and 60, 
respectively yj: 



Table 1. Summary of orbital parameters from Earth fiybys during the last couple of decades, according to Ref. ST}. 



Mission 


Flyby 
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Perigee 


Voo 


A 


Av oof Voo 




Date 




(km) 


(kni/s) 


(mm/s) 


(10-6) 


Galileo 


1990 


2.47 


959.9 


8.949 


3.92 ±0.08 


0.438 


Galileo 


1992 


3.32 


303.1 


8.877 


-4.6 ± 1 


-0.518 


NEAR 


1998 


1.81 


538.8 


6.851 


13.46 ±0.13 


1.96 


Cassini 


1999 


5.8 


1173 


16.01 


-2± 1 


-0,125 


Rosetta 


2005 


1.327 


1954 


3.863 


1.80 ±0.05 


0.466 


MESSENGER 


2005 




2347 


4.056 


0.02 ±0.01 


0.0049 


Rosetta 


2007 




- 2400 




- 




Rosetta 


2009 




2481 




- 





Table 2. List of orders of magnitude of possible error 
sources during Earth fiybys, as discussed in Refs. ^ 
Effect Magnitude 

(m/s^) 



Earth Oblateness 
Other solar system bodies 
Relativity 
Atmospheric drag 
Ocean and Earth tides 
Solar pressure 
Earth infrared 
Spacecraft charge 
Earth albedo 
Solar wind 
Magnetic moment 



10^ 

10-7 
10-7 
10-7 

10-7 
10-7 

10-8 
10-9 

10-9 
10-15 



— K {cos Si — So)- 



(1) 



The constant K is expressed in terms of the Earth's rota- 
tion velocity lue, its radius Re and the speed of light c 
as 



K ^ 



(2) 



suggest that the Earth's rotation may be generating a 
much larger effect than the frame dragging predicted by 
General Relativity. This, however, is in contrast with 
the recent measurements of this effect performed by the 
Gravity Probe B probe [5J, which orbits the Earth at a 
height of about 600 km, well within the onset zone of the 
reported flyby anomaly. 



2. EFFECT ON GNSS SYSTEMS 



In order to discuss the possible use of the available and 
future GNSS constellations to probe this flyby anomaly, 
one should first evaluate to what extent it can affect 
their individual elements. Since the anomalous velocity 
change is only observed before and after fiybys occurring 
at heights of the order of 1000 km, and the GNSS con- 
stellations are in approximately circular Medium Earth 
Orbits (MEO), at about ^ 20000 km, one may empiri- 
cally dismiss any effect. 

One could sharpen the above argument, even though a 
full analysis is impossible due to the lack of spatial res- 
olution and consequent inability to fully characterize the 
spatial dependence of the reported anomaly. A more elab- 
orate discussion is found in Ref. (U. 



This identification is suggestive, as it evokes the general 
form of the outer metric due to a rotating body H, 



with 



Vo 



1 - 2^) idr^ 



(3) 



(4) 



where Vq is the Newtonian potential at the equator. Fol- 
lowing this tentative reasoning, and given the strong lat- 
itude dependence of Eq. ([Til, this expression appears to 



Notwithstanding, one takes as relevant figure of merit the 
anomalous acceleration a ^ IQ-"' m/s^, which may be 
assumed constant for the reasons above. In this case, 
Ref. ||6l shows that no constant acceleration greater than 
10-9 m/s^ can affect the GNSS constellation, since it 
would have otherwise been detected. Thus, one con- 
cludes that the fiyby anomaly, if real, must be due to a 
strongly decaying force, which should drop by four or- 
ders of magnitude with a modest (about than a factor 
four) increase in distance, from r — Re + h 7000 km 
to r ~ 27000 km. As a result, one may safely assume 
that the GNSS constellation is fundamentally unaffected 
by this putative anomaly, and may be thus employed to 
track probes performing gravity assists at the relevant re- 
gion h ^ 1000 km. 
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Figure 1. Velocity error of multi-GNSS tracking of HEO 
spacecraft, adapted from Ref. [12]. Boxes (centered on 
perigee with 4 h width) signal the gap in DNS coverage; 
the horizontal line corresponds to a 20 nim/s accuracy. 



Table 3. Typical accuracies expected from GNSS satellite 
tracking systems for LEO, MEO, Geosynchronous Earth 
Orbit (GEO)and HEO [12. 13. 14. 15]. 





Apogee 


Position 


Velocity 


Orbit 


Height 


Accuracy 


Accuracy 




(km) 


(m) 


(mm/s) 


LEO 


200 to 2000 


10 


10 


MEO 


2000 to GEO 


30 


20 


GEO 


35786 


150 


20 


HEO 


> 35786 


100 


20 



3. GNSS SPACECRAFT TRACKING 



The tracking of spacecraft through GNSS systems is al- 
ready commercially available {e.g. EADS-Astrium's Mo- 
saic Q, NASA PiVoT lis)). These systems are typically 
used to follow satellites in low earth orbit (LEO), at alti- 
tudes below those of the GNSS satellites {h < /iqnss ^ 
20000 km), where the GNSS signal is strongest. 

Nevertheless, the Equator-S mission can receive front 
lobe signal from GPS satellites at an altitude of 61000 km 
|9|. Furthermore, it is worth exploring the possibility of 
using the side and back lobes of the GPS signals ifTOl [TTI 
to establish non-line of sight tracking and avoid the shad- 
ing of the Earth. Clearly, the build up of more constella- 
tions and a rational use of multi-GNSS receivers, able to 
work simultaneously with different systems, will increase 
the accuracy of above-MEO satellite tracking. 

The accuracy of GNSS spacecraft tracking is, under- 
standably, better for lower orbits; however, it should be 
noted that during the apogee of highly elliptical orbits 
(HEO), the velocity is, of course, much slower than close 
to perigee. This allows for the construction of a good 
orbital solution, despite the decreased signal coverage 
lfT2l [T3l . As a result of this trade off, the position and 
velocity accuracies for different types of orbit are some- 
what similar, as depicted in Table |3] 

Recall that there is no full characterization of the anoma- 
lies during the flyby, and these are detected from the mis- 
match between the expected and observed velocities after 
gravitational assist; as stated before, this is due to the in- 
ability of the DSN to track the spacecraft trajectories very 
close to the atmosphere, during a ^ 4 h gap. Regarding 
the possibility of using the GNSS in this region. Fig. [T] 
(adapted from Ref. |12|) shows that, although the veloc- 
ity error is maximum close to perigee, this peak is very 
localized: from a baseline of ^ 20 mm/s during the re- 
maining orbit, it maxes out briefly at ~ 100 mm/s (dur- 
ing the first perigee approach), and converges towards 
^ 50 mm/s in the subsequent perigee passings. By plot- 
ting the aforementioned gap, one sees that accuracies of 
^ 20 mm/s are attainable during approximately half of 
this time interval. 



For the study of the flyby anomaly, one would be inter- 
ested in a high velocity accuracy, at least of the same or- 
der of magnitude as the observed Avoo ^ 1 mm/s. The 
currently available systems provide around 20 mm/s, 
which is clearly insufficient for such study. However, 
the presented accuracies are related to real-time orbit so- 
lutions — which is unnecessary for the purpose of this 
study, and can undoubtedly be improved if offline pro- 
cessing is used, alongside other weak signal tracking 
strategies lfT3l . This, together with the increasing num- 
bers of elements of the available (and upcoming) GNSS, 
lead us to conclude that it is indeed feasible to use the 
latter to test the flyby anomaly — if not presently, then in 
the near future. 



4. PROBING THE FLYBY ANOMALY 



We consider two options to test the flyby anomaly: an 
add-on to an existing mission on a Highly Elliptic Orbit 
(HEO), or a dedicated low-cost mission in either HEO or 
a hyperbolic trajectory. 

In the first option, the choice would be to piggyback a 
GNSS receiver in an existing space mission. Since these 
receivers are relatively small and with reduced power 
consumption 17] H), the host mission could be a small 
low-cost one. At apogee, a highly elliptical trajectory 
would present a comparable (although smaller) velocity 
and height as the reported anomalous gravitational as- 
sists, with the added benefit of allowing for repeated ex- 
periments. 

One can regard as an example the cancelled Inner Magne- 
tosphere Explorer (IMEX) mission of the NASA Univer- 
sity Explorer programme, with a mass of only 160 kg and 
a total budget of 15 M$ (jTe). The IMEX probe was to be 
launched as a secondary payload on a Titan IV launcher, 
but was cancelled due to cost overrun; it would have fol- 
lowed a HEO, as summarized in table |4] which would 
provide a "flyby" velocity at perigee of ~ 10 km/s, close 
to the reported anomalous flybys. 

The more ambitious option of a dedicated mission natu- 
rally has a number of advantages over the former. One 
highlights the choice of orbit that can closely mimic a 



Table 4. Summary of orbital parameters oflMEX's HEO 
and a similar hyperbolic trajectory. 

IMEX Hyperbolic " 



Perigee 
Apogee 
Velocity at perigee 
Eccentricity 
Orbital period 



349 km 
335800 km 
^ 10 km/s 
0.72 

10.5 h 



349 km 

11 km/s 
1.04 



gravity assist — including an hyperbolic one. However, 
as discussed above, a closed orbit of sufficiently high el- 
lipticity would provide for multiple flybys, increasing the 
quality of the obtained data and allowing for a better char- 
acterization of the anomaly. 

Furthermore, possible error sources such as aerody- 
namic and thermal effects close to apogee could be more 
closely controlled with a dedicated mission: for instance, 
the spacecraft could be enclosed in a spherical radio- 
transparent body, so to simplify modelling and reduce 
directional effects; if put into a spin, any accidental 
anisotropics would be averaged out, yielding a much 
cleaner testbed for the desired experiment. 

This mission would require a micro-satellite with a mass 
under 100 kg and a budget caped at less than 15 M$, the 
cost of the IMEX mission. This upper bound is rather 
straightforward to argue by comparison: firstly, one does 
not anticipate any additional spending due to the sim- 
plified spherical design over the more complex IMEX 
probe; secondly, the scientific instrumentation found in 
the latter is replaced by a multi-GNSS receiver in the dis- 
cussed dedicated mission, thus lowering the total cost. 

As the purpose of this paper is to present the feasibil- 
ity of using GNSS to probe the flyby anomaly, this es- 
timative aims at illustrating the low cost of a dedicated 
mission for such purpose. Nevertheless, the actual cost 
could in principle be much lower than 15 M$, not only 
due to the inherently simpler design and instrumentation, 
but also because of the ongoing trend of decreased micro- 
satellite costs — reflecting advances in miniaturization, 
lower power consumption and improved industrial pro- 
cesses ifTTl . 



5. CONCLUSIONS 



When considering the use of Galileo (or GNSS in gen- 
eral) to study the flyby anomaly, one finds that most avail- 
able studies deal with the tracking of spacecraft in real 
time, which is characterized by an insufficient velocity 
and position accuracy to currently detect this hypotheti- 
cal phenomenon. However, since this real time accuracy 
is only one order of magnitude above that required (in 
particular, ^ 10 mm/s vs. ^ 1 mm/s in velocity), one 
expects that this situation could change in the short-term: 



a thorough exploitation of available resources could lead 
to a suitable tracking of spacecraft below the stated ac- 
curacy, by abandoning real time solutions and instead re- 
sorting to offline processing, use of side and back lobe 
tracking, amongst other weak signal tracking strategies. 
Crucially, the use of several GNSS at once should be 
paramount, due to the increased coverage gained from 
the different geometries. 

Thus, one can safely state that there is no a priori issue in 
using GNSS tracking to study the reported flyby anomaly. 
Naturally, this availability is not sufficient, as only space- 
craft equipped with a (multi-)GNSS receiver would allow 
for such study. In this work, we have shown that such 
a mission could be easily deployed, either as an add-on 
package to an existing hub with the required highly ellip- 
tical orbit, or via a dedicated mission. 

While the first scenario would provide a cheap solution, 
we argue that a dedicated mission could be envisaged 
with a higher scientific payoff, while maintaining an over- 
all low-cost approach. 

Regardless of the actual origin of the flyby anomaly (un- 
accounted conventional effect, numerical procedure or, 
more tantalizingly, new physics), we believe it offers a 
low-cost opportunity for displaying the scientific possi- 
bilities opened by the GNSS era — and Galileo in partic- 
ular. 
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